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SUMMARY 


A numerical procedure has been developed to calculate the flow fields 
resulting from the viscous-invlscld Interactions that occur when a strong 
jet exhaust and aircraft flow field coupling exists. The approach used In 
the current procedure is to divide the interaction region Into zones which 
are either predominantly viscous or Invlscld. The flow In the Invlscld 
zone, which surrounds most of the aircraft, is calculated using an 
existing linearized potential flow code. The viscous flow zone, which 
encompasses the jet plume, is modeled using a parabolized Navler-Stokes 
code. The key feature of the present procedure is the coupling of the 

zonal solutions such that sufficient Information is transferred between 
the zones to preserve the effects of the interactions. The zonal 
boundaries overlap with the boundary conditions being the Information link 
between zones. An Iteration scheme Iterates the coupled analysis until 
convergence has been obtained. The procedure has been successfully used 
for several test cases for which the computed results are presented. 
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1.0 INTRODUCTION 


New varieties of military and commercial aircraft of Interest use a strong 
coupling between the jet exhaust flow and the overall airplane flow field 
to enhance performance. These aircraft range from STQL transports such as 
the Boeing YC-14 and the NASA QSRA to combat aircraft with highly 
integrated nozzles that can be used for lift enhancement and thrust 
vectoring. Development of these aircraft requires the ability to achieve 
high performance for low development cost. The complex flow interactions 
characteristic of nozzle Installations with strong flow field couplings 
make the development of such Installations expensive, high risk 
undertakings when using conventional parametric test based design 
approaches. The 3-D flow fields associated with these installations can 
Include mixing layers, wakes, separation, and strong curvature effects as 
illustrated in Fig. 1. With little applicable experimental data available 
and without appropriate analytical tools, the designer of such 

Installations Is faced with a formidable task. 

In the traditional design approach, the designer relies heavily on 
parametric model scale wind tunnel test simulations of the proposed 
installation over the range of nozzle and flight conditions to be 

encountered by the aircraft. The available aerodynamic and propulsion 
performance data base Is used to aid in the selection of a configuration. 
Analytical procedures used are generally 2-0 and, therefore, have limited 
value. High quality model scale wind tunnel tests of complex nozzle 

Installations have proven to be expensive and extremely difficult to 

implement. Often the details of the boundary layers and the jet plume are 
not measured. This leads to uncertain results which cannot be 

meaningfully scaled when designing a larger device. Hence, the test based 
design approach is usually expensive while providing results of unknown 
qua! i ty . 

Developments in computational fluid mechanics and computer technology 
offer the potential for a significant improvement in the design process. 


V. 
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Numerical codes have been developed to model complex 3-0 flows. When the 

appropriate analyses are available, parametric analysis can replace 

parametric wind tunnel testing In the design process. This offers 

distinct advantages since the analytical approach allows the designer to 

maintain precise control over the flow and boundary conditions. Full 

scale Installations can be examined since physical size Is not a 

constraint for a computer program. Furthermore, all of the flow variables 

can be examined In detail to obtain a better- understanding of the flow. 

Wind tunnel testing Is still necessary to validate the analysis and for 

“fine tuning" the design, but the configuration test matrix can be 

substantially reduced In size with an equivalent reduction in development 

• 

cost. As computers Increase In size and speed and codes become more 
efficient, the analytical approach offers an Increasingly practical and 
less expensive alternative for the design process. 

Analysis of the complex flow phenomena associated with the strongly 
coupled nozzle Installations with a single flow analysis would require a 
solution of the complete Navier-Stokes equations. Even with today's flow 
modeling codes and state-of-the-art computers, this is not practical. 
Accurate numerical solutions would require vast amounts of storage and 
days of computer processing time. When the number of parameters to be 
varied in the design process is considered, one can imagine that the 
computational costs would be prohibitive aside from the necessary 
technology developments. 

The problem of numerically modeling these complex flow fields can be 
examined from another perspective. Not all of the flow phenomena of 
interest occur In all regions of the flow domain. One can divide the flow 
domain into separate zones in which only certain phenomena are known to 
exist. The complete Navier-Stokes equations can then be simplified to 
model the flow in each zone. Hence, Instead of one solution of the 
Navier-Stokes equations for the full domain, one has a set of simplified 
zonal solutions that must be carefully coupled together to preserve the 
Interzonal Interactions. This reduction In the overall complexity of the 
problem significantly decreases the computational costs and brings the 
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solution of the problem back Into the realm of what Is practical to 
achieve. The zonal modeling approach to the solution of fluid mechanics 
problems has existed for years. Numerous examples of procedures for 
successfully coupling potential flow codes to boundary-! ayer codes have 
been presented In the literature. These procedures have typically used a 
boundary-layer displacement thickness to couple the effect of the 
developing boundary layer to the potential flow solution. Lemmerman and 
Sonnad (Ref. 1) have recently demonstrated that a surface transpiration 
approach Is equivalent while offering a reduction in computational costs. 
Brune, Rubbert, and Forester (Ref. 2) demonstrated a more sophisticated 
zonal modeling approach by coupling a potential flow code to a 2-0 
Navler-Stokes code to model the flow separation behind an ellipsoid. The 
demonstrated success of zonal modeling for these simpler flows suggests 
that It Is a practical and feasible approach for simulating complex 3-0 
flow fields. 

The objective of the work reported herein was to develop an analytical 
procedure for predicting strongly coupled jet exhaust interactions with 
the overall airplane flow field following a zonal modeling approach. In 
this analysis the flow domain is divided into Invlscid and viscous zones. 
The inviscid zone encompasses the entire aircraft. The viscous zone is 
carved out of the inviscid zone and surrounds the jet exhaust plume which 
Is dominated by viscous interactions. One can envision more zones such as 
the 3-0 boundary layers on the aircraft surfaces or a separation zone on 
the nacelle as depicted in Fig. 1, but these are beyond the scope of the 
present study. In the Inviscid zone, a 3-0 potential flow solution is 
computed by the PANAIR pilot code (Ref. 3). A numerical method for the 
parabolized 3-D navier-Stokes equations (Ref. 4) Is used to calculate the 
flow in the viscous zone. Key features of the present work are the 
development of a solution coupling procedure and an associated iteration 
scheme that result in converged solutions while preserving the necessary 
flow of information between the zones. 
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The following section presents a discussion of the Individual codes. This 
is followed by a discussion of the coupling procedure. Next Is a section 
In which the results of demonstration and validation test cases are 
presented. Finally, the concluding remarks and suggestions for Improving 
the coupled analysis are presented. 


2.0 DISCUSSION OF THE COMPONENT ANALYSES 


To analyze the strong Interactions between the jet plume and the flow 
field surrounding the aircraft, a procedure was developed to couple two 
existing codes. Each of these rades Is particularly suited for predicting 
the flow characteristics In the zone In which It Is used. The 

predominantly Invlscld flow that encompasses the aircraft Is predicted 
using the PANAIR pilot code. The flow In the jet plume Is dominated by 
viscous effects which cause entrainment of the invlscld flow. The flow In 
this zone Is calculated using a code that solves the parabolized 
Navi er- Stokes equations. The analyses Incorporated in these codes are 
discussed In the following sections. 

2.1 PANAIR PILOT CODE 

The PANAIR pilot code (henceforth called PANAIR) Is a preliminary version 
of an advanced panel code intended to solve a variety of boundary value 
problems In steady subsonic or supersonic invlscid flow. This code Is 
discussed In Reference 3. Oust the subsonic capability will be considered 
here, since the coupled analysis procedure is at present limited to 
subsonic flow. 

The PANAIR analysis is based on the assumption that for a wide range of 
flow conditions, a potential flow solution will substantially describe the 
flow past a prescribed configuration. This Implies that the flow is 
assumed to be inviscid and irrotatlonal . 

The flow past the configuration Is considered to be a small perturbation 
of a uniform flow that exists far upstream. PANAIR generates the 
potential flow solution In terms of the perturbation velocity potential, 
<&, which Is used to calculate the local perturbations in the flow 
properties. The basic boundary condition employed In PANAIR is that the 
total mass flux vector be parallel to solid surfaces. Non-solid surfaces 
can have mass fluxes normal to the surface. 


The method used in PANAIR to compute the potential flow is to superimpose 
fundamental solutions of the partial differential equation for 4. These 
fundamental solutions are sources and doublets whose locations are 
prescribed but whose strength must be determined to meet the prescribed 
boundary conditions. The sources and doublets are distributed on 
continuous networks of quadrilateral “panels" that approximate the surface 
of the configuration, as well as the vortex sheets shed from trailing 
edges and other surfaces such as inlet faces and jet plumes. Each panel 
Is divided Into subpanels on which the source and doublet strengths are 
approximated by linear and quadratic functions respectively. These 
functions contain the unknown source and doublet strengths which are then 
determined by the simultaneous solution of the algebraic system of 
equations. The resulting potential and velocity fields are then 
determined. The pressure field can then be calculated, and the forces and 
moments on the configuration .'an be computed. 

2.1.1 Panel ing 

Although a detailed discussion of paneling concepts and problems is 
presented In Reference 3, some of the salient features are repeated here 
to emphasize the need for understanding the correct method for paneling a 
given configuration for use In the coupled analysis. 

In general a configuration is divided Into a system of networks which are 
composed of panels. The number of networks depends on the complexity of 
the configuration. The number of panels In a given network will depend 
upon the detail and accuracy desired from the analysis. To use PANAIR as 
part of the coupled analysis, networks for the surface of the 
configuration and Its wakes must be supplied. Networks defining the jet 
plume and its wake are required. Additional networks are necessary for 
obtaining the flow properties at points off the body. These panels which 
do not disturb the flow, are specified with zero jumps In the potential 
and the normal component of mass flux. 

The network surfaces are defined by an array of grid points which provide 
the coordinates for the corners of the quadrilateral panels. The array Is 
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arranged In rows and columns. Care must be exercised In organizing these 
rows and columns since they. In effect, define the upper and lower 
surfaces for each network. Suitable diagnostics are provided by PANAIR to 
aide In debugging the network arrays. When networks abut, It Is necessary 
that common grid points appear In all of the abutting network arrays. If 
the abutting networks have similar grid densities, the grid points used 
along the abuttment must be exactly the same In all of the arrays, unless 
one or more of the networks Is free of doublets. This will Insure that 
the doublet strength will be continuous at the network boundaries. 
Doublet networks must also abut along complete edges such that their 
network corner points coincide. This requirement often necessitates the 
division of a configuration Into numerous small networks when junctures, 
such as between a wing and strut, are encountered. Failure to abide by 
the above rules will lead to disastrous results. 

2.1.2 Boundary Conditions 

The panels are covered with a contl nuors distribution of sources and 
doublets. Boundary conditions are Imposed at selected "control points" on 
the networks to generate the necessary set of equations which are solved 
for the unknown singularity strengths. Networks can have just sources or 
just doublets in which case the boundary conditions for the other 
singularity type are not required. Special options exist for the 
treatment of wakes. 

The user determines the flow behavior by specifying the types of boundary 
conditions on each of the networks. A great deal of flexibility exists, 
and Reference 3 describes the available options. The most commonly used 
types of boundary conditions are those for impermeable surfaces, the 
Mori no-type boundary conditions. These are used on the wetted surfaces of 
an aircraft. For the coupled analysis the jet plume network has boundary 
conditions that allow the specification of the perturbation mass fluxes 
normal to the panel surfaces. Thus, the plume network Is a permeable 
surface that entrains fluid from the surrounding potential flow. The use 
of this boundary condition is discussed in detail in Section 3.0. 
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2.1.3 Inputs and Outputs 

The inputs to PANAIR are divided Into four basic categories: 

(1) General specifications 

(2) Flow conditions 

(3) Network geometry specifications 

(4) Boundary condition specifications. 

All of the Inputs with the exception of the geometry are generally 
specified through FORTRAN coding in a special input subroutine. The 
geometry and part of the plume boundary conditions are provided In Input 
files for the coupled analysis. 

The outputs can be controlled to a large degree by the user. Options 
exist for printing out various types of diagnostics that allow a thorough 
examination of the health of the Inputs before the solution routines are 
actually executed. Local aerodynamic data is always printed for every 
center-control point on every network. The output subroutine Is modified 
for use In the coupled analysis to generate extra data files which 
t,rm<-’er boundary conditions to the viscous zone. 

2.2 PARABOLIZED NAVIER- STOKES CODE 

The parabolized Navier-Stokes code ( PMS) was developed to calculate a 
particular class of three-dimensional compressible viscous flows. The 
flows of interest are characterized by a predominant flow direction. In 
this class of flow, downstream disturbances have a negligible influence on 
the upstream flow. Thus, the assumption is made that the propagation of 
perturbations In the upstream direction by convection, diffusion, or 
pressure can be neglected for this particular class of flows. Cross 
stream propagation Is allowed and in fact can be significant. Flows of 
this type are generally classified as being "parabolic". They exist in 
numerous internal and external flows ranging from complex 3-0 duct flows 
to 3-D jets mixing with a freestream. 
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, The time averaged Navler-Stokes equations with suitable closure for 

turbulent flows are acknowledged as being sufficient for analyzing complex 
viscous flow fields. When the above parabolic approximations are 
Implemented, these equations are simplified so that they can be solved in 
a practical design oriented computer program. The Navler-Stokes equations 
are parabolized by neglecting the streamwlse diffusion terms and 
decoupling the streamwlse and cross stream pressure gradients. The 
streamwlse pressure gradient Is assumed to be uniform at each cross- 
sectional station. The resulting equations are elliptic In the cross 
stream and are often considered parabolic In tne streamwlse direction, 
however, this Is more historical than formal. In the PNS code a simple 
marching. solution procedure was implemented which eliminates communication 
with the downstream flow while solving the parabolized Navler-Stokes 
equations. Other features ~*f the PNS code, which are discussed In tne 
following sections, are the use of a general coordinate system, the 
transformation of the flow equations, the pressure-continuity relations, 
the turbulence model, the solution procedure, and the features that were 
added to automate the coupling with PAHAIR. 

/ 

2.2.1 General Coordinate System 

Flows In geometrically complex domains are difficult to compute on 
standard orthogonal computational meshes such as Cartesian or cylindrical 
coordinate systems. The difficulties arise because the mesh points do not 
naturally fall on the boundaries of the flow domain of Interest. This 
results In cumbersome differentiation and Interpolation schemes at the 
boundaries that can achieve only low levels of accuracy. For many flows 
it Is desirable to maintain high accuracy at the boundaries. For this 
reason significant effort in recent years (Refs. 5-8) has been aimed at 
generating general coordinate systems that are fitted to the boundaries of 
the flow domain. In the PNS code a curvilinear mesh is formed in which 
the boundary mesh points always fall on the natural boundaries of the flow 
domain and the Interior mesh points conform to the boundary shape. This 
allows the use of standard differencing expressions and maintains high 
levels of accuracy at the boundaries. The boundary- fitted mesh is 
particularly suited for coupling with a panel method which also can handle 
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arbitrary geometries. Numerous methods exist for automatically generating 
the boundary-fitted meshes. These methods offer varying degrees of 
complexity and capability. Orthogonal meshes can be generated, but 
absolute orthogonality has not been demonstrated to b«s necessary or worth 
the effort, though highly skewed meshes are undesirable. 

The boundary-fitted £, n » a coordinate system used In the PNS code 
can be' generated using one of two options. The first method generates the 
interior mesh by simply interpolating between the Inner and outer 
boundaries. This method has the capability of concentrating the mesh 
points near the outer boundary. The second method Is similar to that 

discussed by Thompson, et al (Ref. 5). This method requires the solution 

of a coupled set of elliptic quasl-llnear partial differential equations. 
A significant amount of flexibility exists since the Interior mesh can be 
controlled by the use of suitable forcing functions. This ultimately 
allows one to automatically Increase the mesh density in regions where 
greater accuracy Is required. The 3-0 mesh Is generated in a step-by-step 
fashion. As the flow solution Is marched from plane to plane the mesh at 
the next downstream plane is computed as the immediate precursor to the 
solution of the flow equations. Thus, the computational mesh for the 

complete flow domain only exists when the flow Is completely predicted. 

2.2.2 Transformed Flow Equations 

To use the general boundary-fitted computational mesh, the flow equations 
are formulated with 5 , n, and <7 as the independent variables. Using 
standard transformation relations the Cartesian primitive variable form of 
the Navler-Stokes equations is transformed to the 5 , n » and a 
coordinate system. These equations are parabolized with the a 

coordinate specified to be the streamwlse or predominant flow direction. 
The transformed equations are further simplified by requiring that the 
cross-sectional planes be parallel and perpendicular to the z coordinate. 
This limits the amount that the mean centerline of the flow domain can 

deflect from a straight line, but any significant curvature would 
typically indicate a non-parabolic flow. In 5 , q, and a coordinates 
the parabolized steady, three-dimensional continuity, momentum and energy 
equations are as follows: 
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f The appropriate diffusion terms, which are too cumbersome to present, are 

| 

| represented by the S on the rlghthand side of each equation. The 

Cartesian velocity components have been retained. The T which 
I appears in Eq. (Id) represents the constant streamwise pressure gradient 

j for a given cross-section. In addition to the field equations, an 

i equation of state is included, and a calorically perfect fluid is 

assumed. Sutherland's formula is used to relate the viscosity, u, to 

r r- 

l 

| temperature. 

i 2.2.3 Pressure-Continuity Relations 

Although the streamwise and lateral pressure gradients have been decoupled 
j to parabolize the flow equations, they are implicitly coupled through the 

* continuity equation. The functional purpose of the pressure gradients is 

j to insure the conservation of mass both globally and locally. Equations 

have been developed to guarantee this effect, since mass conservation is 
an essential requirement for the numerical simulation. A relation for 

r is derived from the w-momentum equation and an integral definition 

f for the mass flow rate in the diffuser. Mass is globally conserved by 

adjusting P" and Iterating the w-mcmentum equation until the mass flow 
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rate reaches a specified level of accuracy. Global mass conservation is 
required before local conservation can be considered. 

Local continuity Is satisfied indirectly by means of the lateral pressure 

gradients P„ and P . These pressure gradients are calculated such 
5 n 

that the u and v velocity components generated by Eqs. (lb) and (lc) will 
simultaneously satisfy the continuity equation (la) everywhere In the 
computational domain to a specified tolerance. The method used In the PNS 
code Is similar to methods developed for time-dependent flows (Refs. 9,10). 

Starting with a cell A In the lateral plane with vertices 1, 2, 3 and 4, 
an Iteration scheme Is developed. The u- and v -momentum equations 
Indicate that a perturbation of the cell pressure has the following effect 
on the surrounding velocities. 


U 1 ■ “1 - W^A 

(2a) 

v i “ v i ■ t $ jjw 

(2b) 

-i. 

a 

1 



where the ~ represents the previous values. Equations (25 are substi- 
tuted Into the discrete form of the continuity equation (la) to yield. 


AP A 3 uf (3) 

where C is the discrete continuity equation for cell A, 8 is a combination 
of mesh derivatives and velocity coefficients which result from the above 
substitution, and w is an over-relaxation parameter i Using (3) aP a 
is calculated. The surrounding velocities are then updated using 
equations (2) and the pressure P A is updated. This procedure Is 
repeated at each computational cell. The computational plane Is iterated 
until C decreases below a specified small value at all cells. A 
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relaxation factor of 1.5 seems to be optimal as long as the C distribution 
Is relatively smooth; otherwise, slight under-relaxation Is necessary 
until the smoothness develops. This pressure iteration scheme does not 
directly Include the momentum convection and diffusion effects, however 
iteration with the momentum equations does implicitly produce these 
effects. Since iteration is already part of the solution procedure, no 
additional Iteration is required. 

2.2.4 T urbulence Model 

Closure of the set of flow equations for turbulent flows is achieved by 
modeling the Reynolds stresses that appear in the time-averaged 
Navier-Stokes equations. In the PNS code these stresses are modeled by 
replacing the laminar viscosity p with a turbulent eddy viscosity 
p . The eddy viscosity is calculated from the turbulence kinetic 
energy k and the turbulence kinetic energy dissipation rate e. 

\i - C pk /e 
t u 

The values for k and e are calculated using two additional transport 
equations (Ref. 11), which have the same form as Eq. (le) with additional 
terms that model production and dissipation. 

2.2.5 Boundary Conditions 

The three velocity components u, v, and w are set to zero on walls. An 
adiabatic wall boundary condition Is used with the energy equation. To 
minimize the mesh requirements in the lower regions of boundary layers, 
law-of-the-wall functions are used to represent the velocity distribution 
near the wall. These functions are used to calculate the wall shear 
stress and the production and dissipation source terms in the turbulence 
model equations for the wall region. The use of wall functions allows k 
and e to be accurately predicted in the vicinity of a wall. 

2.2.6 Numerical Solution Procedure 

All of the differential equations presented above have been transformed to 
finite-difference form by standard second order centered difference 
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approximations for 5 and n derivatives and by one-sided upstream 
differences for a derivatives. Since this makes the numerical 

approximation formally second order for a suitable streamwise step size,, 
numerical noise can always be expected to generate spurious solutions that 
can be catastrophic unless It Is properly controlled. A noise filter 
based on the work of Forester (Ref. 12) has been Incorporated in the 
numerical algorithm to reduce the effects of the computational noise to 
negligible levels. 

An iterative marching solution procedure was developed to allow converged 
numerical solutions after just one pass through the flow domain. An 
Initial set of data is required at the starting plane. A solution Is 

obtained at each successive cross section before a step Is taken 
downstream to the next lateral plane. The solution of the nonlinear flow 
equations at each plane is achieved by an Iterative ADI procedure. 
Iteration Is required at each plane to reduce the accumulation of 
truncation errors that can result from linearization. Iteration allows a 
fully converged implicit solution at each plane. Usually, just a tew 
iterations are required to achieve reasonable accuracy unless the flow is 
changing rapidly In the streamwise direction. The flow chart in Fig. 2 
illustrates the basic solution procedure used in the PNS code. 

2.2.7 Features Added for PANAIR Coupling 

To facilitate a fully automatic iterative coupling procedure for the PNS 
code and PANAIR, It was necessary to add several capabilities to the PNS 
code. A routine was written to generate the panel networks that define 
the jet plume boundary and the boundaries of the viscous flow domain. The 
mass fluxes through the boundaries are calculated in an additional 
routine. A special procedure was developed to track the jet plume with 
the computational mesh such that the mesh could be used efficiently. The 
vortlcity of the plume is monitored to locate the plume edge. With this 
known, the position of the mesh at the next state can be projected. This 
also ensures that the panels defining the plume can be positioned outside 
of the mixing region where viscous effects dominate. 


2.2.8 Inputs and Outputs 

The Inputs to the PNS code can*be divided Into four categories: 

(1) Controls and options for mesh generation 

(2) Physical constants 

(3) Boundary points for the mesh 

(4) Initial conditions for the flow variables 

The inputs are all read from a data file. The flow variables can be 
initialized internally for a few special cases, but they are generally 
read from the data file. 

The main output to the program provides the calculated values for each of 
the flow variables at all mesh points. Diagnostic outputs are available 
to provide a more detailed history of the evolution of the flow. An 
output file can also be generated for use in computerized graphic displays. 
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3.0 COUPLING PROCEDURE 


The interactions between the jet plume and the surrounding flow field can 
be correctly simulated by the proper coupling of the component analyses. 
It Is not sufficient to run each code once in its separate zone. The 
codes must interact with each other In such a manner that the solution for 
the entire flow field Is convergent and unique within the limits allowed 
by the algebra Incorporated In the codes. Each code must provide 
Information to the other code that adequately describes the physical 
processes being modeled in that zone. The PNS code, which calculates the 
jet plume development, must provide PANAIR with the effects of entrainment 
at the plume boundary. PANAIR must use this information and the other 
aerodynamic effects of the aircraft to calculate the potential flow field 
and in turn provide boundary conditions for the PNS code. The procedure 
developed for coupling PANAIR and the PNS code is discussed in the 
following sections. 

3.1 OVERLAPPING COMPUTATIONAL ZONES 

The domain in which the flow field is to be calculated is divided into 
computational zones for the Individual codes. The Inviscid zone Is 

oriented to include all regions which can be adequately modeled by the 
PANAIR potential flow solution. The viscous zone is positioned to 

surround the jet plume where the viscous interactions are predicted by the 
PNS code. The location of the boundaries for these zones is a significant 
aspect of the coupling procedure. 

Abutting the zones such that the boundaries are coincident does not 
provide the necessary flow of information for a convergent coupled 
analysis. Since both codes in effect solve boundary value problems, a 

specification of boundary conditions on the coincident boundaries would 
lead to unique solutions in each zone that are functions of those boundary 
conditions. Unless a priori knowledge of the flow properties on the 

boundaries was available, the flows predicted in each zone would not 
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necessarily bear any resemblance to the physical flow. New Information 
that could be used to update boundary conditions would not be available 
since the solutions are boundary condition dependent and their boundary 
conditions are Identical. Iteration Is useless without new Information. 
Hence, coincident boundaries lead to a coupled analysis that does not 
allow the codes to Interact and exchange Information and yields an 
Iteration procedure that will not be reliably convergent. 

These considerations led to the conclusion that for the proper coupling of 
these two codes, the boundaries of the computational zones eannot be 
coincident. Furthermore, the boundaries of each zone should be arranged 
such that the boundary conditions are dependent upon the solution in the 
neighboring zone. This allows the necessary transfer of Information from 
one zone to another. In the present coupling procedure, this Is 
accomplished by overlapping the zonal boundaries. This is Illustrated in 
Fig. 3 for the simple case of an axlsymmetrlc nozzle exhausting Into a 
co-flowing freestream. Note that the boundaries of the viscous zone 
extend well past the edge of the jet plume and into the region that Is 
inviscld flow. The inviscid zone boundaries extend Into the viscous zone 
to a close proximity of the plume edge where the flow is still dominated 
by inviscld effects. The overlapping region is shared mutually by both 
zones. Therefore, the boundary conditions applied in the PANAIR code on 
the plume boundary can be derived from the flow properties calculated by 
the PNS code. Similarly, the boundary conditions for the viscous zone can 
be obtained from the PANAIR solution. The overlapping boundaries 
therefore provide the necessary communi cation link between the two zones 
that allows a converged coupled solution for the flow domain. 

The solutions computed by PANAIR and the PNS code are consistent within 
their numerical accuracies for the overlapping region. If the flows under 
consideration are restricted to those that are parabolic, the P^ term 
in the Equation (Id) will be negligible for external flows. By realizing 
that the flow In the overlapping region Is effectively inviscid, the flow 
equations reduce to $he Euler equations. This region is also Irrotational 
since no mechanism Is present to generate vorticity. In the overlapping 
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region the PNS code Is actually solving the potential flow equation which 
is the fundamental equation of the PANAIR analysis. Therefore, the two 
solutions should be consistent In this region. Due to the various 
approximations Incorporated Into these two codes, exact consistency cannot 
be achieved, but the accuracy will be comparable to that obtained for the 
complete coupled solution. The overlapping region should be kept 
reasonably small, since It would be inefficient to have a large region In 
which the flow Is in effect calculated twice. 

3.2 BOUNDARY CONDITIONS 

i 

The specification of the boundary conditions on the overlapping boundaries 
Is straightforward. The solution from the PNS code is used to calculate 
the perturbation mass flux through each of the panels that form the jet 
plume boundary network. This is accomplished by averaging the predicted 
velocity components and density on the corners of each panel and using the 
result to calculate the component of the mass flux vector that Is normal 
to the panel. The array of perturbation mass fluxes is stored for access 
by the PANAIR code. 

The boundary conditions for the PNS code are determined from the PANAIR 
solution In a similar manner. Special panel networks that do not disturb 
the flow are positioned in the Invlscld zone to form the boundaries for 
the viscous zone. The predicted velocities on the panels are used to 
generate the boundary conditions for the PNS code. In the PNS code, these 
velocities are also used to compute the mass flux entering the viscous 
zone. The upstream boundary of the viscous zone, Fig. 3, forms the 
initial data plane for the PNS code. On that boundary, the panel 
velocities are used to calculate initial conditions, therefore, that 
region Is initially Irrotatlonal . 

3.3 MESH AND PANEL GENERATIQJ 

The generation of the panels on the overlapping boundaries Is directly 
related to the generation of the computational mesh for the PNS code. All 
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other panels are directly Input Into PANAIR through the normal Input 
procedure. Various options can be used to define the boundaries for the 
overlap region. Each of these options affects the method used for 
generating the mesh and panels. To facilitate the discussion of these 
options, consider the diagram in Fig. 3. 

The outer boundary of the overlap region (and of the viscous zone) Is 
designated to have a computational mesh index of J 3 N. The Inner boundary 
of the overlap region (the plume boundary for the Inviscld zone) has the 
index J=»JNP. In the PNS analysis, the JNP mesh points are computationally 
like all other Interior mesh points In the viscous zone. The correct 
usage of the coupled analysis requires that the JNP boundary lies beyond 
the edge of the viscous plume while not crossing the J a N boundary. 

The simplest option for generating the N and JNP mesh points is to specify 
the N boundary points at the initial plane and project these points to 
each downstream plane such that all planes in the viscous iuiie will have 
identical boundary contours. The JNP points are specified at the initial 
plane on the edge of the plume and projected downstream with a given slope 
such thait the JNP boundary expands. The panels on the N and JNP 
boundaries are constructed by using the mesh points on these boundaries as 
the comer points for the panels. The initial plane panels are similarly 
generated from the initial plane mesh points. This option has limited 
capabilities since the expanding JNP boundary could cross the N boundary 
at some point which would cause the PNS code to blow up. This problem is 
rectified by allowing the N boundary to expand with the same slope as the 
PNS boundary. The boundary slope must be large enough for JNP boundary to 
capture the plume and small enough to efficiently use the available 
computational mesh. Only simple plumes can be considered since the plume 
trajectory is limited to a narrow path. 

The complexity of the shapes of the plumes to be calculated is enhanced by 
another option in which the N boundary mesh points are specified for the 
Initial plane and for all succeeding downstream planes. The JNP points 
are specified in a like manner. The associated panels are generated as 


discussed above. Clearly, this option requires that the user have some 
knowledge of the expected trajectory of the plume. This Information is 
often available when the analysis Is being used to parametrically examine 
a case for which some experimental data exists for similar flows. Care 
must be used In sped flying the JNP boundary to guarantee It surrounds the 
jet plume. 

Another option was developed for the coupled analysis to automatically 
track the jet plume. At the Initial plane the JNP boundary points are 
specified. The N boundary points are computed based on the JNP points and 
the number of mesh points one desires to have between JNP and N. At each 
of the downstream planes the JNP points are positioned using the vortlclty 
In the most recent plane to locate the plume edge. This procedure 
Involves calculating the vortlclty at all points and determining where the 
vortlclty approaches zero, which Indicates the edge of the plume. This 
data Is used to project the JNP points dcwnstfv&fg to positions that will 
be outside of the plume. The maximum change In position can be controlled 
by specifying a maximum slope from upstream to downstream. The N boundary 
points are then calculated as before, and the panels are generated. This 
approach to generating the mesh Is flexible since the boundaries of the 
viscous zone can move and distort with the plume. The viscous zone can 
then be compact and computationally efficient. 

3.4 ITERATION PROCEDURE 

One Iteration of the solutions In the invlscld and viscous zones will 
generally not be sufficient to yield a converged solution for the total 
flow dnmain. The essential ingredient Is the information transfer between 
the two zones. This Is accomplished by repetitively updating the boundary 
conditions and generating new solutions. An iteration procedure was 
developed for the coupled analyses to automatically perform the necessary 
data manipulations and transfers without requiring Interruptions by the 


user. 


To start the Iteration process, an Initialization of the boundary 

conditions on the overlapping boundary of one of the zones Is required. 
This Is done by making an “educated'' guess as to what the boundary values 
should be. The better the guess the fewer the number of Iterations to 
reach convergence. The zone In which one chooses to Initiate the 

iteration is not Important, unless auxiliary Information Is available to 
yield a better estimate of the boundary conditions In one of the zones. 
If the viscous zone Is the starting point, the velocity components on the 
Initial plane and the outer boundaries must be provided. The freestream 
velocity components are often adequate for this Initialization purpose. 
By starting with the Invlscld zone*, one Is faced with the task of 
estimating the perturbation mass fluxes on the plume p/inel network. First 
the overlapping boundary networks must be constructed by executing the 

panel generation portion of the PNS code. Unless knowledge of the jet 

entrainment Is available, the Initial plume boundary fluxes can only be 
guessed. Using a constant flux Is sufficient to start the Iteration 
procedure. 

As a demonstration of the iteration procedure, assume that the solution Is 
started In the Invlscid zone and refer to Fig. 4. The overlapping 
boundary panels are generated by the PNS code. Then the perturbation mass 
fluxes on the plume panel network are Initialized. PANAIR is executed In 
the invlscid zone. The potential flow solution is used to calculate the 
boundary conditions for the viscous zone. The next step Is the execution 
of the PNS code In the viscous zone. This will Include the generation of 
new overlapping boundaries if the plume tracking capability Is utilized. 
New values for the boundary conditions on the plume panel network are 
computed from the viscous zone solution. At this point these new boundary 
values can be compared with the old values. If a prescribed level of 
convergence Is achieved, the iteration procedure can be terminated. 
Otherwise, the iteration loop will continue for a given number of 
iterations. The alternative of starting the solution In the viscous zone 
has the same Iteration loop with the exception that the starting point 
would be the fourth step in Fig. 4 In which an Initialization of the 
viscous zone boundary conditions Is made. 


4.0 VALIDATION CASES 


To demonstrate and validate the usefulness of the coupled PNS and PANAIR 
codes, several test cases were set up and run. The main purpose of these 
tests was to demonstrate the coupling procedure and show that converged 
solutions can be achieved. In addition Information about the behavior of 
the coupling procedure as the solutions Iterate In the different zones 
could be obtained and used to understand and Improve the existing 

procedure. While the accuracy of the overall solution depends on the 
accuracy obtainable by each code, the accuracy of the coupled codes can be 
validated to some level of satisfaction by comparing the computed 

solutions with existing experimental data when possible. The following 
cases provide a variety of the applications of the coupled analysis. 

4.1 AX I SYMMETRIC JET IN A FREESTREAM 

The first test case, a subsonic axisymmetrlc jet in a freestream, was used 

primarily as an aid in the development of the coupling procedure. The 

simplicity of this case, see Fig. 3, was desirable since the effects of 
the coupling procedure and the overlapping boundary condition could be 
readily isolated. In addition experimental data and previous 
computational work were available for comparison. To minimize the 
computational costs, the symmetry or the case was used to limit the 
solution domain to half of the jet. In addition just five mesh points 
were used on the cl rcumference of the jet, while 21 points were usee 
radially. In the overlapping region seven mesh cells were used to 
separate the Inner Invlscid zone boundary from the outer viscous zone 
boundary. The plume tracking option was used to locate the plume edge by 
monitoring the vortlclty. This guaranteed that the number of mesh cells 
In the jet plume would be the same at all stations since the mesh expands 
at the same rate as the plume. It also ensured that the panels generated 
for the plume network would be outside of the plume where the strong 
viscous Interactions dominate. 
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The freestream for this case had a Mach number of M»0.2. The ratio of the 
jet core velocity to the freestream velocity was 2.17. In the jet, the 
initial velocity was given a constant value which gave the Initial 

velocity profile In the viscous zone a top hat distribution. This type of 

distribution has large velocity gradients which typically give rise to 

computational noise. The noise filtering Incorporated in the PNS code 
effectively controlled the noise In the Initial regions of the viscous 

zone until the steep velocity profiles had diminished. The turbulence 
quantities were initialized with values that were empirically determined 
for developing free shear layers. A constant total temperature was used 
In the jet and freestream. The Initial marching step size, aZ, for the 
PNS code was equal to 0.25 of the initial jet radius. The step size was 
then gradually Increased to a maximum of 1.25 Initial jet diameters for 
one run and 2.5 diameters for another to examine the effect of step size. 
The coupled analysis was started by generating the panel networks In the 
overlapping region, initializing the perturbation mass flux on the plume 
network, and executing PANAIR in the invlscld zone. Then the zonal 
solutions and the overlapping boundary conditions were iterated using the 

9 

automated coupling procedure. 

Exparlmental data for the axlsymmetric jet is available for the decay of 
the jet centerline velocity and the spreading rate of the jet. Rodl (Ref. 
13) has also computed axisymmetrlc jets using a turbulence model like that 
incorporated In the PNS code. 

The results computed for this case using the coupled PANAIR/PNS analysis 
are presented in Figs. 5 and 6. The centerline velocity decay Is 
underpredicted. Fig. 5, when compared with the experimental data, but the 
spreading rate calculations. Fig. 6, compare more favorably. However, the 
present calculations are similar to Rodl's. Rodi concluded that the 
two-equation turbulence model was adequate for calculating the 
axisymmetrlc jet when the velocity ratio was large. However, he also 
found that the accuracy of the calculations is sensitive to initial 
conditions, mesh density, and step size. The sensitivity to aZ is 
confirmed by the present results. 
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Of greater Importance was the Information obtained about the coupling 
procedure. The plume tracking option Is effective In minimizing the 
required amount of computational mesh by placing the mesh and panels in 
the proper locations. However, since the final plume boundary Is not 

known until the coupled solution converges, the plume tracker, which 
calculates the new boundary as the solution Is marched through the viscous 
zone, moves the boundary with each Iteration of the PNS code. In effect 
the overlapping boundaries approach their final positions as the flow 
solution converges. This has a detrimental effect which tends to reduce 
the overall convergence rate. The boundary conditions for the viscous 
zone are calculated at the previous boundary positions, however, they are 
used In the PNS code with the new boundaries. The movement of the 
boundaries has associated area changes and normal velocity errors which 
lead to an error in the mass flux Into the viscous zone. It is this small 
error In entrainment that tends to reduce the convergence rate. This 
effect Is demonstrated In Fig. 7 which shows the convergence histories for 
the boundary conditions at several locations along the developing jet. At 
the upstream locations the convergence rate Is substantially greater than 
at the downstream locations. When the plume tracking option was 
activated, a maximum allowable spreading rate for the plume edge was 
specified. In the initial region this maximum slope was achieved for 
several steps during each Iteration through the viscous zone. Hence, the 
boundary did not move with Iteration or affect the local convergence 
rate. Farther downstream the maximum slope was not achieved, the boundary 
moved during each Iteration, and the convergence rate was reduced. 

Convergence of the viscous zone outer boundary position is not required 
for convergence of the flow solution, but the plume boundary needs to be 
converged. To reduce the detrimental effect of a moving outer boundary, 
the standard plume tracking option could be modified to move this boundary 
for a couple of iterations until a rough position is found. Then the 
outer boundary should be frozen to allow a more rapid convergence of the 
flow solution. 
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A related point of interest is the effect of convergence on the streamwise 

pressure gradient term F. As discussed in section 2.2.3, F, Is 

adjusted to satisfy overall continuity. For the axl symmetric jet, F 

should always be zero. However, errors in the boundary conditions lead to 

an inconsistent mass flow rate which results in a nonzero F. As the 

__ a 

solution converges these errors diminish and approaches zero at all 

streamwise stations. 

4.2 AX I SYMMETRIC JET OVER AN AIRFOIL 

The main purpose for developing the present coupled analysis was to 

provide an analytical method for predicing the effect of the jet plume on 
the aerodynamics of a given aircraft configuration. In section 4.1 the 
coupled PANAIR/PNS analysis was demonstrated for a simple jet to provtde 
convergent solutions with an accuracy that depends on the limitations of 
the component codes. Jets of substantially greater complexity could be 
examined at various angles of attack, but that would not provide 
additional insight into the jet/aircraft interaction problem. To test the 
coupled analysis for the interaction problem, a simple test case was set 
up for an axi symmetric jet positioned over an airfoil of constant cross 
section. This case provides a degree of flexibility since the jet 
position can be varied to determine whether the analysis predicts the 
expected trends. 

A schematic diagram for this test case is depicted in Fig. 8. The 

axisymmetric nozzle of diameter, D n , is positioned a height, h, above 
the airfoil. The leading edge of the airfoil and the nozzle exit plane 
coincide. The half span of the airfoil was 2.5D n> The initial 
conditions for the jet and the freestream conditions were identical to 
those discussed above for the isolated axisymmetric jet. The bilateral 
symmetry plane was used to simplify this problem. The plume tracking 
option was also used for this case. A fairly coarse computational mesh 
and panel network representation was used to conserve computer costs since 
only trends were to be predicted and no experimental data was available. 
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The coupled procedure was initiated In the Invlscid zone with a PANAIR 
solution and then iterated. This case was run with the nozzle at three 
different heights above the airfoil. 

The basic result obtained from this study was the net lift on the 
airfoil. The airfoil drag was not considered since the boundary layers 
were not part of this investigation and the jet was not allowed to scrub 
the wing. Fig. 9 presents a plot of the airfoil lift versus the height of 
the nozzle. The lift has been normalized by the lift of the airfoil 
without a jet present, and the nozzle height is normalized by the nozzle 
diameter. One can observe that the increase in lift is greatest when the 
nozzle is closest to the airfoil and the interaction is the strongest. 
The entrainment of the jet increases the circulation about the airfoil 
which enhances lift. Fig. 9 can only be used to observe trends since the 
actual increment in lift will be a function of the length of the airfoil 
among other things. However, the trends predicted by the coupled analysis 
agree with what is found experimentally. This offers encouragement for 
the continued use of this analysis. 

The effect of the jet/airfoil interaction on the convergence of the 
coupled flow solution was substantial. The convergence history for the 
boundary conditions on the plume network panels closest to the airfoil 
surface is presented in Fig. 10, with h/D =0.9 the perturbation mass 
fluxes approach their final values smoothly and quickly in just three 
iterations. The convergence behavior was substantially different with 

h/D =0.775. The rate of convergence was still reasonable, but four 

iterations were necessary, and it appears that one more iteration may have 
been warranted. Part of this slowdown is attributed to the movement of 
the jet plume boundaries since the plume actually was deflected upwards 
and then downwards as it progressed over the airfoil. The stronger 
interaction at the lower height also affects convergence since the plume 
has a greater influence on the PANAIR solution for the airfoil which in 
turn has a large influence on the PNS solution for the plume. In general 
large perturbations reduce the rate of convergence, particularly when the 
plume axis must move to a converged position. 
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4.3 USB NOZZLE/WING 


A case of practical interest is that of a rectangular nozzle blowing over 
the upper surface of a wing, see Fig. 11. This upper surface blowing 
(USB) type of nozzle has been successfully used in recent years on such 
aircraft as the Boeing C-14 and the NASA QSRA. In the present case the 
computed results are compared with the experimental investigation 
discussed in Reference 14. 

As shown in Fig. 11 the nozzle exit plane Is positioned approximately 0.32 
of the wing chord, C, from the leading edge. The nozzle half width, R n> 
is .21C. The half span of the wing is 6.4 R n „ Since the coordinates 
for the wing cross section were not available, they were estimated from a 
diagram in Reference 14. The angle of attack of the wing was set to zero; 
however, for the experimental data the wing had an angle of attack which 
is not known. Therefore, only qualitative comparisons with the data were 
attempted. 

The flow properties were set up to simulate one of the experimental test 
conditions as closely as possible. The freestream Mach number was 
M »0.6. The freestream total pressure and total temperature were 1400 
lb/ft^ and 520°R. The velocity profile at the nozzle exit plane was 
calculated from the experimental total pressure surveys. The average 
total pressure in the nozzle was greater than the freestream. Part of the 
nozzle flow was slightly supersonic which is beyond the capability of the 
PNS code, however, most of the viscous zone was subsonic which allowed the 
PNS code to successfully generate a solution. The boundaries of the 
inviscid and viscous zones were fixed since the plume tracking option was 
not designed to work on a solid surface. The computational mesh in the 
viscous zone was too coarse to provide more than qualitative results. 

The principle result from this test case is the effect of the nozzle on 
the pressure distribution on the upper surface of the wing. Since the 
wing was analyzed at zero angle of attack and the wing cross section was 
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just an estimate, the calculated pressure distribution did not resemble 
the experimental data, and a comparison would be of little valde. The 
predicted trends on the other hand were Identical to those found 
experimentally. Figure 12 provides the calculated chordwlse pressure 
distribution at a position Y=1.8 R n from the nozzle centerline (PANAIR 
coordinates are used). The distribution for the wing alone and the 
wing/nozzle with the jet are shown. It Is observed that with the jet the 
pressure drops below the wing alone curve at approximately X/OQ.2 and 
rises above it at X/C=0.6. Hence, the effect of the jet plume is to 
reduce the lift of the wing at the spanwise position between X/C of 0.2 
and 0.6 and increase the lift beyond X/C of 0.6. This trend Is exhibited 
In the experimental data Including the points at which the curves cross 
(Fig. 16 In Reference 14). This favorable comparison offers encouragement 
for the potential of the coupled analysis. Since part of the nozzle flow 
Is supersonic, the pressure distribution along the nozzle centerline at 
Y/R n a 0 cannot be accurately calculated. By modifying the method for 
calculating pressure in the PNS code, these supersonic regions could be 
predicted. In addition this case should be run again with denser 
computational meshes and panel networks to assess the effects of 
truncation errors. 

4.4 623 V/STQL MODEL 

The Grunman 623 V/STOL wind tunnel model has been extensively tested by 
NASA/ Ames (Reference 16). The 623 model Is a fighter- type configuration 
Fig. 13, with the nacelle Integrated Into the wing at the wing/body 
junction. The nozzle Is an ADEN rectangular cross-section nozzle. Fig. 

14, with the thrust deflector flap blended into the wing trailing edge. 
The inboard edge of this flap forms a juncture with the body when the 
deflector angle is 0°. The flap can be moved to deflect the exhaust 
plume to increase lift at low speeds and provide increased maneuverability 
for the vehicle. At zero deflection angle the exhaust plume positioned at 
the wing trailing edge can still provide Increased lift at low freestream 
Mach numbers since the jet entrainment induces additional circulation 
about the wing. 


29 


The 623 model has been tested for free stream Mach numbers ranging from 
0.4 to 0.9 and angles of attack of nominally 0° to 13°. Cruise and 
combat variations of the ADEN nozzle were used. The thrust deflector 
angle ranged from 0° to 20°. The experimental data consists primarily 
of pressure coefficient distributions on the wings, nacelles, and body and 
of force and moment measurements . Exhaust plume measurements were not 
obtained. Limited total pressure data were taken upstream of the nozzle 
throat, and a few static pressures were measured along the nozzle walls. 
The nozzle stagnation conditions were set to simulate: (1) no flow, (2) a 
flow through nacelle, and (3) a supersonic exhaust plume. Since the PNS 
code does not currently have the capability for calculating supersonic 
plumes, a limited number of the 623 model test cases can be of uss for 
demonstrating the coupled analysis. However, It was felt that analyzing a 
complete aircraft configuration would be beneficial for understanding the 
coupling procedure. In addition, the 623 model had already been 
successfully paneled by NASA Ames. The paneled 623 model Is shown In 
Figure IS. The sparseness of the panels was required for the case to run 
on the Ames CDC 7600 computer. The additional panels needed by the 
coupled analysis necessitated the use of the Boeing CRAY computer. 

4.4.1 Internal Nozzle Flow Analysis 

To obtain a complete set of flow conditions at the nozzle exit plane for 
use In the coupled analysis, the internal flow field for the nozzle was 
computed. By choosing the appropriate options in the PNS code, the 
boundary conditions are switched to model an impermeable, no-slip, 
adiabatic wall. In this mode the program has been used to successfully 
calculate numerous internal flows (Reference 4). 

The 623 model ADEN nozzle has two basic throat configurations, cruise and 
combat. The combat configuration, which was used for the present 
calculations, has a larger throat area. The internal geometry of the 
nozzle. Fig. 16, starts with a circular cross-section downstream of a 
choke plate assembly. This Is followed by a rapid transition to a 
rectangular cross-section. At this point the nozzle width is fixed, but 
the height varies through the throat section to the exit plane. 
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The Initial flow conditions downstream of the choke plate are required as 
Input to the PNS code In order to calculate the internal flow field for 
the nozzle. The only available measurements from the test data were the 
mass flow rate and the total temperature and total pressure, hence, 
certain assumptions about the Initial flow field were made. Since the 
choke plate tends to remove Irregularities In the flow, uniform profiles 
for velocity, temperature, and the turbulence quantities were assumed to 
exist at the Initial plane. The average velocity and static pressure were 
calculated from the known flow conditions. Using this set of initial 
conditions and the adiabatic, no-slip wall boundary conditions, the. PNS 
code was used to march the solution to the nozzle exit plane. At this 
point a file containing the calculated flow field data was saved to be 
used for the plume Initial conditions when the coupled analysis was used 
for the external flow field. At pre^nt it Is not possible to couple the 
Internal flow to the external flow, since a fully elliptic viscous 
analysis would be required. 

4,4.2 External Flow Analysis 

The paneling for the 623 model with the thrust deflector at zero degrees 
deflection angle was available from NASA Ames. Since the panels for other 
deflector angles were not available, the coupled analysis could not be run 
to demonstrate the effects of thrust deflection. However, the zero angle 
data could still be used for comparison with the numerical results. 

4. 4.2.1 623 Model with No Nozzle Flow 

In preparation for Implementing the coupled analysis the paneled 623 model 
was analyzed using PANAIR alone. The experimental nozzle pressure ratio 
NPR, was 1.0 for this case. Indicating that no air was flowing through the 
exhaust nozzles. This was modeled by covering the nozzle exit with a 
panel network having an Impermeable boundary condition. The freestream 
Mach number was 0.6 and the angle of attack was zero degrees. The purpose 
of this baseline case was to determine the adequacy of the PANAIR solution 
in the regions of Interest without the interactions caused by the jet 
plume. The PANAIR results were compared with the experimental data at the 
stations where y a 0.11, 0.41, and 0.76 of the aircraft's semi -span 
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distance, b/2. Here, y Is the PANAIR spanwlse coordinate direction. 
These stations correspond to locations on the upper surfaces of the 
fuselage, nacelle and wing. For this and all other 623 model cases the 
wake panel networks were angled up at five degrees. 

The pressure coefficient, C p , distribution along the upper fuselage 
station Is provided In Fig. 17. Clearly, the PANAIR solutions follows the 
experimental data reasonably well. The most significant discrepancies 
occur on the aft 20 percent of the fuselage. In this region the boundary 
layer could be thick making separation a possibility, which would account 
for the substantial differences. The boundary layer at this station was 
not investigated experimentally. The other discrepancies may be 
associated with an Insufficient number of panels to provide adequate 
resolution or with bad data points. 

The results for the nacelle upper surface are shown In Fig. 18. The 
discrepancies between the PANAIR solution and the experimental data at 
this station are significant only at locations between 60 and 80 percent 
of the nacelle length, where the predictions seem to be out of phase, and 
at the leading edge. These differences are probably related to a lack of 
panels and viscous effects. 

The wing Cp distributions are compared in Fig. 19. The agreement 
between the analysis and experimental data Is good except in the vicinity 
of the leading edge. This again may be due to an inadequate number of 
panels to resolve this region of strong acceleration. 

Overall, the PANAIR solution for this non-flowing nozzle case seems to 
provide reasonably accurate results. Similar results were obtained at 
NASA-Ames using the same paneling. By increasing the panel density, 
better agreement with the experimental data may be ootained, but that is 
beyond the scope of this Investigation. 
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4. 4. 2. 2 623 Model With Flow Through Nacelle 

The 623 model was tested at nozzle total pressures approximately equal to 
the freestream total pressure. This had the effect of simulating a flow 
through nacelle since the nozzle flow total temperature was equal to that 
In the freestream. The NPR for this case was 1.22. The freestream Mach 
number and angle of attack were 0.6 and 0°, respectively. Although the 
effect of the plume Is small, the test data Indicates It is noticeable, 
Fig. 20. Before this case could be analyzed using the coupled PANAIR and 
PNS codes, several modifications to the paneled 623 model were necessary. 

Since the plume was exiting the nozzle at a low pressure ratio. It would 
not expand. The analysis from Section 4.4.1 was used to calculate the 
exit velocity profile, which had a mean value near the freestream velocity 
as expected. The lack of a significant velocity ratio between the 
freestream and the plume suggests that virtually no mixing would take 
place on the deflector ramp. In view of this the ramp was modified by 
enclosing It with panels that had the effect of moving the nozzle exit te 
the deflector trailing edge. This greatly simplified the complexity of 
the coupling precedure. 

A further modification to the nozzle was required in order to accommodate 
the overlapping boundaries of the coupled analysis. Since the edge of the 
deflector flap In the 0° position has a juncture with the fuselage, part 
of the overlapping mesh needed by the PNS code at the Initial plane would 
be positioned inside the aircraft. This unallowable situation was 
remedied by adding a short impermeable extension to the nozzle. Since the 
aft end of the fuselage rapidly converges, a one Inch extension was 
sufficient for the mesh boundaries to avoid the body. Though the 
extension does modify the plume flow slightly, the negligible expansion 
and mixing In this region make the effect insignificant. 

The initial plume flow conditions for the PNS code were calculated as 
discussed In Section 4.4.1. The initial velocity profile had a peak value 
that was below the local external flow field. Therefore, this simulated 
flow through nacelle plume becomes a wake that interacts with the aircraft 
flow field. 
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The panels that define the nozzle trailing edge extension were modeled 
using just sources as opposed to the source/ doublet panels used on the 
aircraft body. This was done to avoid problems associated with the 
requirement for a continuous doublet distribution at the Intersection of 
the trailing edge wake panels with the deflector flap. The panels 
generated by the PNS code to define the plume were also specified to be 
pure source panels. 

The coupled analysis was Initiated by using the PNS code to generate the 
plume shape and associated panels. A uniform perturbation mass flux was 
specified on the plume for the first Iteration of PANAIR. Each code was 
Iterated four times, which was sufficient for convergence. Each Iteration 
of the PNS code resulted in a new set of plume panels and mass flux 
distribution. After each Iteration of both codes the results were 
examined before the next Iteration was allowed to proceed. Convergence of 
the overall solution for this case was monitored by plotting the Cp 
distributions for each Iteration at the fuselage, nacelle, and wing 
stations. These distributions seemed to rapidly converge. 

The results for this case are presented in Figure 21-23. The predicted 
Cp distributions show the same trends as the experimental data. The 
Incremental change in C p Is greater for the coupled analysis solution 
along the aft half of the fuselage and nacelle. The C p distribution 
along the wing shows a favorable comparison between the analytical and 
experimental data. The agreement is particularly good along the aft end 
of the wing. The overall solution appears to be modeling the effect of 
the plume on the aerodynamics of the aircraft reasonably well. 

4. 4. 2. 3 623 Model with Jet Plume 

Test data for cases with a jet plume were taken at nozzle pressure ratios 
that make the jet supersonic. Since the PNS code was not designed for 
supersonic flow, a case was set up for a plume with the Initial Mach 
number set to 0.95. The intent of this case was to determine whether a 
jet could have a strong effect. A uniform velocity distribution at the 
exit plane was used. The ratio of the jet velocity to the freestream was 
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1.55. The total temperature and static pressure were the same as for the 
previous case. Four iterations of the coupled analysis were run for this 
case. 

The predicted results for this case are shown In Figures 24-26. Clearly, 
the Increase in C p Is minimal. This Is not totally unexpected. The low 
velocity ratio does not lead to significant mixing of the plume with the 
surrounding flow, which means that the entrainment effects should be 
small. At the relatively high freestreara Mach number of this case, the 
effects of entrainment cannot be felt upstream as they are for low speed 
flows. An examination of the test data for the high NPR cases reveals the 
same trend demonstrated with the analysis. The small Increase in 
with increasing NPR demonstrates that the coupled analysis procedure 
allows the jet entrainment effects to be calculated by the PNS code and 
transmitted to the PANAIR code through the overlapping boundaries. Test 
data with larger influences from jet entrainment on the wing circulation 
would be preferable for better analysis validation. 
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5.0 FUTURE WORK 


As a result of the work discussed In this report and the need for advanced 
numerical procedures, several recommendations for future work are 
presented In this section. These recommendations are aimed at increasing 
the generality and efficiency of the present coupled analysis. Extending 
the analysis to allow the prediction of more complex flows can be pursued 
within this framework as more powerful analyses become available. 

5.1 IMPROVEMENTS TO THE PRESENT COUPLING PROCEDURE 

The present coupled procedure Is still expensive for flow field solutions 
which require a substantial number of iterations for convergence. Methods 
for Improving the convergence rate need to be Investigated. The plume 
tracking option needs further development to include a provision for 
freezing the outer viscous zone boundary after the general plume location 
is established. This would stop the Introduction of mass flux errors at 
the boundary that Impede convergence. Panel and mesh density studies can 
then be accomplished jsore efficiently. 

5.2 SUPERSONIC VISCOUS JET ANALYSIS 

As discussed In Section 4.3, the PNS code will not accurately predict 
supersonic flows. However, nozzles with supercritical pressure ratios are 
often encountered, and it is desirable to have the capability to predict 
the supersonic jet plumes. The PNS code can be modified to implement 
pressure relations that are valid for supersonic flows. This would have 
to be done such that the resultant analysis is still compatible with the 
overall coupling procedure. 
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6.0 CONCLUSIONS 


The most significant conclusion of the work discussed in this report is 
that a successful procedure has been developed for coupling PANAIR and 
parabolized Navler-Stokes solutions to allow numerical simulations of the 
strong interactions between jet plumes and the overall aerodyamlcs of 
aircraft configurations. Overlapping Invlscid and viscous computational 
zones provide the means for Information transfer between the solutions. 
The coupling procedure is automated and yields convergent solutions. This 
zonal procedure was applied to four test cases; an axi symmetric jet in a 
co-flowing stream; an axi symmetric jet above an airfoil; a rectangular jet 
on an airfoil; and a V/STOL fighter airplane configuration with 2-0 
nozzles at the wing trailing edge. Good quantitative agreement with 
experiment data was obtained for the first case. For the remaining cases, 
qn^’l -j ■ha + 'jya ggrggmgn* wac Qnmnirfc|Ki TSSUl tS tQ data Indicate 

that the panel densities for the PANAIR solutions and the mesh densities 
for the Navler-Stokes solutions significantly affect the predictions. 
Additional mesh refinement studies heed to be performed to evaluate the 
procedure and to establish panel and mesh density criteria. 

Both of the component analyses are limited as to the types of flows that 
can be simulated. The present coupling procedure can be extended by 
modifying the existing analyses and adding new analyses. The complexity 
of the flows to be predicted will determine which extensions need to the 
made. The technology developed for the present procedure will be directly 
applicable to future work. 
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